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Abstract —The mitigation of nonlinear distortion caused by 
power amplifiers (PA) in Orthogonal Frequency Division Multi¬ 
plexing (OFDM) systems is an essential issue to enable energy 
efficient operation. In this work we proposed a new algorithm 
for receiver-based clipping estimation in OFDM systems that 
combines the Iterative Hard Thresholding method introduced 
in f5| with the weighting corresponding to the estimated proba¬ 
bility of clipping used in (Sj. Thereby a more general amplifier 
input-output characteristic is considered, which is assumed to 
be unknown at the receiver side. Further, we avoid the use 
of dedicated subcarriers and formulate the recovery problem 
solely on reliably detected sub-carriers. Through simulations, we 
show that the proposed technique achieves a better complexity- 
performance tradeoff compared to existing methods. 

I. Introduction 

With the recent leap from 3G to 4G, the amount of power 
consumption in mobile devices has increased considerably. 
The 4G operates on Long Term Evolution (LTE) in which 
Orthogonal Erequency Division Multiplexing (OEDM) is used 
for the wireless communications. However, the use of multi¬ 
carrier modulation such as OEDM leads to the higher Peak-to- 
Average Power (PAPR) and thus higher power consumption in 
an RE Power Amplifier (PA) as well as the Digital-to-Analog 
(ADC) converters in the mobile devices. In a typical 4G 
system, RE PAs have to satisfy stringent linearity requirements 
and thus accounts for a significant amount of the total power 
consumption. Therefore, in order to operate the amplifiers 
at maximal possible efficiency and reduce the requirements 
on the DAC, it is important to allow a partial clipping of 
the OEDM signal, while envisaging at the same time the 
compensation of the resulting distortion at the receiver side. In 
fact, when operating at such lower input back-off for the PA 
results also into a certain out-off-band radiation and a certain 
in-band distortion. This paper consider the compensation of 
the in-band distortion. 

In this context, the work Q investigated the iterative 
maximum likelihood (ML) detection of the clipped OEDM 
signal based on the turbo decoding approach and developed 
an computationally efficient algorithm called power amplifier 
nonlinearity cancellation (PANG). The main disadvantage of 
this method is that the input-output characteristic of the PA 
has to be perfectly known at the receiver side. 

In anther recent line of work, a sparse reconstruction al¬ 
gorithm, called Support Agnostic Bayesian Matching Pursuit 
(SABMP) has been proposed in IS), IT], [|4], |8| for the 


recovery of clipped OEDM signals has been investigated based 
on the sparse nature of the clipping error while choosing a 
certain number reliable sub-carrier for Setting up the sparse 
recovery problem. Nevertheless, it has been observed, that 
this algorithm has relatively higher complexity. Therefore, 
a new sparse reconstruction approach based on Weighted 
Iterative Hard Thresholding (WIHT) is proposed which has 
significant lower computation complexity than that of SABMP 
for almost the same performance. Moreover, a more general 
amplifier characteristic is introduced and analyzed, which is 
assumed to unknown at the receiver. Einally, the proposed 
sparse reconstruction algorithm is applied and simulated for 
a single-input single-output (SISO) system. 

Our paper is organized as follows. Section |II] describes 
the general system model including the amplifier model 
and formulates the recovery problem. In Section |III] the 
proposed Weighted Iterative Hard Thresholding algorithm is 
presented. Section |IV] compares the computational complexity 
of this technique with two other widely used methods, namely 
Support Agnostic Bayesian Matching Pursuit (SABMP) 
and Power Amplifier Nonlinearity Cancellation (PANC) El- 
In Section |V] we provide some simulation results to reveal 
the usefulness of the method, and Section |VT] gives a brief 
conclusion and suggests some areas for future work. 

II. System Model 

Let us consider an OEDM system in which the incoming 
information bits are mapped onto L-ary QAM constellation 
and concatenated to form an A-dimensional frequency domain 
data symbols vector X G The time-domain vector x 

can be obtained by x = F^X, where F is a A x A DET 
matrix ||T|. 

The time-domain signal x has high PAPR due to the 
addition of different frequency sub-carriers. Because of this, 
the power amplifier (PA) will operate in the non-linear region, 
introducing clipping for the higher values of x. The clipping 
model we are considering is a generalization of the soft 
clipping used in ID in that it also introduces phase distortion 
in the non-linear region. We motivate this by the fact that the 
behavior in this region is generally unpredictable and cannot 
always be compensated, while phase distortion in the linear 
region is usually mild and can be compensated with pre¬ 
distortion. The samples of the clipped signal Xp are then given 


by: 


\^j{a,rsx{i)+4>{\x(i)\)) ’ 

where t is the clipping threshold and (j){\x{i)\) is a phase 
distortion function defined as the argument of the Band 1 PA 
model given by 3GPP in lO. We note that our proposed sparse 
recovery technique is able to recover the clipping without any 
prior knowledge of this phase distortion function. 

The clipping can be seen as the addition of a sparse signal 
c in the time-domain signal x and can be written as 

Xp = X -b c, (2) 

which can be written as follows 

Xp = + c. (3) 

In order to avoid the ISI at the receiver, a cyclic prefix (CP) 
is appended to the time-domain signal. This CP is removed 
at the receiver and then FFT is performed at the receiver to 
obtain the received signal as 

y = Hxp -I- z, (4) 

where z is the AWGN noise with the variance cr^. H is the 
circulant channel matrix due the CP insertion/removal property 
and is given as H = F*^AF. The received time-domain signal 
y is then written as follows 


y = F^AFxp -b z. (5) 

The received signal in the frequency domain can then be 
written as 

y = A{X + C) + Z, (6) 

the estimated received signal is obtained by equalizing the 
channel as follows 

x = x-^y = {x + c) + x-^z, (7) 

In order to determine the unknown vector c, X needs to 
be eliminated. The process for elimination is carried out by 
projecting onto a set of reliable carriers. This is the set of 
carriers where we have more certainty that the distortion 
introduced by noise and clipping was not big enough to move 
the symbol to a different decision region in the constellation. 
The reliable carriers are chosen according to a reliability 
measure that depends on both the distance to the closest 
constellation point and the angle to it: 



Fig. 1. Reliability of observation relevant to distance and angle to different 
constellation points 


was taken from a, and is motivated by the fact that a 
deviation from the constellation point in an angle close to 
G {0,1, 2,3} makes the subcarrier less reliable than a 
deviation in an angle close to {2k + l)f, k G (0,1, 2, 3}, due 
to the fact that the distance to the second closest constellation 
point is higher in the second case. In summary, samples that 
have a shorter distance to their closest constellation point, or 
whose angle to it is closer to {2k -b 1)^, k G (0,1, 2, 3} will 
get a higher reliability value. The P subcarriers with highest 
value of TZ{X — {X)) are chosen as reliable. 

Taking the difference between the estimate X and the 
associated vector {X) after the sheer and then projecting onto 
the chosen reliable carriers yields 


J{X-{X))= J{{X + C- {X))) + JA-^Z (9) 

= JC + JA^^,Z (10) 

= JFc-bJA~^,Z, (11) 

where J is an P x N binary selection matrix indicating the 
locations of the reliable carriers. Above equation can then be 
written as 


n{X - {^)) = jX {X)\ ^ 


y/2d 

min 

y/2d 

min 


cos(46'_^_^_^^-b7r), (8) 


where dmin is the distance between two constellation neigh¬ 
bors, {X) is the closest constellation point to X, and 
is the angle to it, as depicted in Fig. [T] This geometric measure 


y = Ac + Z', (12) 

where y = J{X - {X)), A = JF and ,2' = JA-^Z. 

This is a compressed sensing model that can now be 
processed for the unknown c by using any CS algorithms. 
In the following, an efficient CS algorithm will be used 
for determining the unknown sparse vector c. Note that, in 
contrast to m, we do not include the phase information into 







the sensing matrix A and recover a complex sparse vector c 
in order to be able to compensate also for phase distortion in 
the non-linear region of the PA. 

III. Weighted Iterative Hard Thresholding 
(WIHT) 

The IHT algorithm IS) is an efficient greedy technique for 
the sparse recovery problem. However, it cannot be directly 
applied for the CS model in (O and some modifications 
are required in order to make use of this robust method. We 
call our proposed modified algorithm Weighted Iterative Hard 
Thresholding (WIHT). 

The summary of the proposed technique is shown in Al¬ 
gorithm [T] The algorithms takes X, y, A and the estimated 
number K of active taps in c as inputs. We note that the 
performance of our method does not degrade noticeably if 
this value is overestimated, therefore, as proposed in 1^, we 
choose K to be slightly higher than the number of elements 
of |A^3^| that are greater or equal than half its maximum 
value. 


K = 


j : \afy\ > i ||A^y|| 


(13) 


where aj denotes the j-th column of A. The clipping 
threshold r is approximated by taking the maximum absolute 
value of the estimated clipped signal Xp = F^X. Afterwards, 
a weighting vector w is calculated to approximate the proba¬ 
bility that each sample was clipped. As this probability must 
be higher for samples closer to the clipping threshold, and 
be constrained to the range [0,1], a negative exponential of 
the difference between the approximated threshold f and the 
estimated clipped signal value was chosen as the weighting 
function ||8]: 


w = exp{-(f-|xp|)}. (14) 

W is the diagonal matrix with the weighting entries along 
its diagonal. The function Hk{-) as shown in the algorithm 
below can be defined as 


Hk{A, b) 


A^b, A^b > max/f (A^b) 
0 otherwise, 


(15) 


where maxfc(-) gives the fc-th maximum value. In other words, 
the function iJ/f(A, b) is used to find out the k maximum 
values of the correlation between the matrices A and the vector 
b. The remaining correlation values which are not included in 
the k maximum values are then forced to zero. In the algorithm 
below, the correlation is computed between the measurement 
matrices A and the observation vector y and then the support 
of the sparse vector is computed by simply taking the support 
of the output of the function 

Up until now, a single iteration of standard IHT has been 
run. We found out that the values of the taps of the sparse 
vector c are generally not well approximated after only one 
iteration, but the support of usually coincides with that 


of c. Therefore, Weighted IHT takes the support iSiht of 
and performs a Best Linear Unbiased Estimate (BLUE) of c 
over it; 

g = (16) 

where As,^^ = is obtained by taking the columns of 

A whose indices are in the support 5 iht- This BLUE estimate 
approximates the conditional expectation over the observa¬ 
tion and the computed support, c Ri E [c|y, iSiht] ■ Einally, 
the actual pre-clipped transmitted signal can be obtained by 
subtracting the estimated sparse vector c from the estimated 
clipped signal Xp. 


Algorithm 1 Weighted Iterative Hard Thresholding (WIHT) 

Input: X,y, A, K 
Estimate of clipped signal: Xp = 

Threshold estimation: T = lIXpIloo 
Weighting vector: w = exp{— (r — |xp|)} 

Weighting matrix: W = diag{w} 

Support of Sparse Vector: 5 iht = Supp(iJif (WA^y)) 

BLUE Estimate: c = E[c|y,5iHT] ~ 

Output: X = Xp — c 


IV. Computational Complexity 

The computational complexity of the reconstruction algo¬ 
rithms compared in this paper was estimated as the number of 
complex multiplications they need. 

The complexity of the proposed Weighted IHT algorithm 
comes mainly from the calculation of the pseudo-inverse 
(AgAs) AgV, which is O (p . 

The W-SABMP algorithm needs to compute this pseudo¬ 
inverse several times, but the efficient implementation pro¬ 
posed in Section IV of m reduces its complexity to O {P |5|) 
by exploiting the previous results. Even then, the calculation 
needs to be done N times for each possible support size 
up to |5|, making the overall complexity of this algorithm 
O 

Einally, the PANC algorithm performs an PET and an lEET 
of size N in each iteration. We found that two iterations 
are needed to achieve BER results comparable to the other 
two algorithms, and therefore PANC has a complexity of 

0{N\ogN). 


V. Simulation Results 

The proposed Weighted IHT algorithm was tested and 
compared with the existing techniques W-SABMP and PANC 
in terms of both uncoded BER and computational complexity. 
The performance without recovery and with the oracle-LS 
method (where the support of the clip signal c is perfectly 
known at the receiver and a least-squares solution is applied 
to estimate its values) are also given for reference. Several 







experiments were carried out to evaluate the effect of the varia¬ 
tion of different scenario parameters on these two performance 
measures. 

All the experiments simulate an OFDM system with N = 
512 subcarriers and 16-QAM modulation. The clipping ratio 
is defined as: 

CR=—, (17) 

where is the standard deviation of the input signal to the 
PA, and t is the clipping level. The clipped OFDM signal goes 
through a 4-tap random complex channel (generated using a 
Gaussian distribution), and is affected by AWGN. The SNR 
is defined at the receiver input: 


SNR = 10 log 
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(18) 


All the clipping estimation algorithms are applied after channel 
equalization. The selection of reliable carriers is made ac¬ 
cording to the geometric reliability measure introduced in 121, 
which is computationally faster than the Bayesian approach 
(0{N) instead of 0{NM), where M is the QAM modulation 
order), while getting almost equally good results. 

We used the efficient implementation of the W-SABMP 
algorithm proposed in Section IV of fh). The PANG algo¬ 
rithm 13 is run for two iterations, as one was not enough 
to achieve results comparable to the other algorithms. All 
algorithms estimate the clipping level as the absolute value 
of the maximum sample of the estimated received signal in 
the time domain. 
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Fig. 2. BER vs P for CR = 1.3 and Ei,/Nq = 15 dB. 



A. Experiment I: Performance vs number of reliable carriers 

For the first experiment, the clipping ratio was fixed at 
CR = 1.3, and the Ei,/Nq at 15 dB, and the performance of 
the algorithms was evaluated for varying number of reliable 
carriers P. 

Fig. 12] shows the results. As expected, the execution time 
of Weighted IHT and W-SABMP increases linearly with P. 
Our proposed technique executes in the order of N times 
faster than the efficient implementation of W-SABMP, and 
has a comparable speed to that of PANC. In terms of BER, 
the sparse recovery algorithms benefit from an increase of the 
number of reliable subcarriers, up to a point where too many 
subcarriers with errors are incorrectly selected as reliable and 
the performance degrades again. The best results seem to be 
around P = 275, so we chose this value for the remaining 
experiments. 

B. Experiment II: Performance vi Eb /iVg 

With a clipping ratio of CR = 1.3, and P = 275 reliable 
carriers, the performance of the three algorithms was measured 
for different values of E^/Nq, obtaining the plots in Fig. [3| 
The BER results for Weighted IHT are only slightly worse 
than those for W-SABMP, with a computational complexity 
two orders of magnitude lower. This confirms that the pro¬ 
posed weighting pi = exp{— (7 — ip [i])} greatly increases 
the convergence speed of these two algorithms, making the 


Fig. 3. BER vs Eij/Nq for CR = 1.3 and P = 275. 

benefits of the other factor p (3^|iS) marginal, and in most 
cases not worth the considerable additional computational cost 
at which they come. 

Our technique also obtains comparable or better BER results 
than PANC, especially at high SNR, where the reliability 
measure is more likely to be correct. The speed of both 
algorithms is almost the same with these parameters, and will 
be better for Weighted IHT if N is reduced (because |iS| will 
decrease proportionally) or if the clipping ratio increases. Even 
with unfavorable conditions (high N and low clipping ratio), 
the difference in speed is very small, and tolerable given the 
fact that Weighted IHT does not need to have the PA model 
of the transmitter, which is not readily available in the uplink 
of a mobile communications channel. 

C. Experiment III: Performance vs CR 

In this experiment, we aimed to check the range of values of 
the clipping ratio for which our proposed algorithm gives good 
results. We fixed the Eb/N^ at 15 dB and P = 275 reliable 
carriers, and plotted the BER and runtime over CR, as shown 
in Eig. m Again, W-SABMP performs only marginally better 
than Weighted IHT in terms of BER, supporting our idea that 
the weighting is the main contribution to their convergence 
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Fig. 4. BER vs CR for E},/No = 15 dB and P = 275. 


rate. The results are also comparable to PANC, slightly out¬ 
performing it at severe clipping levels and somewhat worse at 
higher CR. 

As for the execution time, a strong inverse dependence with 
the clipping ratio can be observed for the sparse recovery 
algorithms, due to the quadratic relationship between their 
computational complexity and the size of the support set 
of the clip signal |5|, which decreases when CR increases. 
Weighted IHT is still two orders of magnitude faster than W- 
SABMP, and from the two graphs we can verify that it is also 
advantageous with respect to PANC when the clipping ratio 
is low (no higher than 1.4). 

VI. Conclusion 

An weighted iterative hard thresholding algorithm for mit¬ 
igating the nonlinear effects of power amplifiers is proposed. 
Compared to the existing PANC technique 121, our algorithm 
obtains similar results, with slightly better BER performance 
and slower execution for milder clipping, and vice versa for 
more severe clipping. However, unlike PANC, our technique 
is able to recover even phase distortion without knowing the 
PA model, which is not available in the uplink of a mobile 
communication cell. In this scenario. Weighted IHT with sim¬ 
ilar performance but fewer assumptions can be advantageous 
with respect to the existing algorithms. 

The proposed technique also leaves some areas open for 
future work, which are mainly the optimization of the weight¬ 
ing function and of the number of reliable carriers to select 
depending on the parameters of the scenario. 
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